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Abstract 

We examine the effects of higher derivative corrections on rj/s, the ratio of shear viscosity 
to entropy density, in the case of a finite i?-charge chemical potential. In particular, we work 
in the framework of five-dimensional N = 2 gauged supergravity, and include terms up to four 
derivatives, representing the super symmetric completion of the Chern-Simons term ylATr (RAR). 
The addition of the four-derivative terms yields a correction which is a effect, and in general 
gives rise to a violation of the rj/s bound. Furthermore, we find that, once the bound is violated, 
turning on the chemical potential only leads to an even larger violation of the bound. 
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I. INTRODUCTION 

Over the past decade the development of the AdS/CFT correspondence has 
led to a new way of thinking about strongly coupled gauge theories. Although the original 
and best studied example of the AdS/CFT duality connects Af = 4 supersymmetric Yang- 
Mills to type IIB string theory on AdSs x 5*^, the duality has been extended to a variety of 
cases, and can describe confining gauge theories with features that are qualitatively similar 
to QCD. In recent years the AdS/CFT correspondence has proven to be a valuable tool 
for better understanding thermal and hydrodynamic properties of field theories at strong 
coupling. In particular, it has been applied to the realm of heavy ion collisions, with the 
aim of providing a more realistic description of the strongly coupled quark-gluon plasma 
(QGP). 

In the context of RHIC physics, a quantity that has played a special role is the ratio of 
shear viscosity to entropy density, rj/s (see e.g. ^ and references therein). Weak coupling 
calculations in thermal field theory predict t]/s ^ 1, while elliptic flow measurements 
at RHIC seem to indicate a very small ratio < ry/s < 0.3, showing that the QGP 
behaves like a nearly ideal fluid, and is in the strong coupling regime. Motivated by such 
observations, there has been a large effort to apply AdS / CFT methods to the calculation 
of various transport coefficients. The AdS/CFT "program" is particularly valuable given 
that lattice methods (which work well for equilibrium, or thermodynamic, quantities) fail 
for non-equilibrium processes. 

Furthermore, developments resulting from the AdS^FT correspondence prompted 
Kovtun, Son and Starinets (KSS) to postulate a bound ^] for rj/s, according to which all 
fluids would obey 

s - in ^ ^ 



The bound, which seems to be satisfied by all substances in nature, was later shown 
to be saturated in all gauge theories with a dual supergravity description in the large 
and A = Qym-^ limit. Moreover, the universal value rj/s = l/47r ~ 0.08 falls into the 



2 



experimental range observed at RHIC. Finite A corrections to the leading supergravity 
result were explored in Q], which considered curvature terms of the form ~ a'^R^ in 
Type IIB supergravity on AdS^^ x . The result was that the leading finite A corrections 
increase the ratio in the direction consistent with the bound: 



7] 1 
s An . 



l + 15C(3)A-=^/2 _ (2) 



However, r]/s bound violations were subsequently observed in the presence of curvature 
squared terms @, I^, Q, ll|. In the context of the AdS/CFT correspondence, such terms 



correspond to finite N corrections and lead to 



!?=ifi-£^i (3) 



s 47r \ a 

where a and c are the central charges of the dual CFT. Thus, violation will occur provided 
c — a > 0. The central charges are known to be equal in the large limit [13], with 
a = c = 0{N'^), but differ for finite A^. For the supergravity examples studied so far, 
the leading corrections on the CFT side lead to c — a > 0, implying violation of the 
bound by finite N effects [12]. (It is an interesting question on its own to ask whether 
one can have string theory constructions whose dual description allows for c — a < 0.) 

In this paper we investigate what happens to the rj/s ratio in the presence of non-zero 
chemical potential. In particular, we focus on the chemical potential corresponding to 
turning on a U{1)r background of the M = 2 system. To leading order in the supergravity 
approximation, the i?-charge chemical potential does not affect the calculation of ril s, as 

n n n 

was shown in |14J . 1151 . Ilq . However, it is interesting to examine whether this is still the 
case once higher derivative corrections are included. Furthermore, if rj/s is affected by 
i?-charge, it would be useful to see whether the KSS bound violation gets larger or smaller 
function of chemical potential . 



^ Ideally, one could imagine tuning the chemical potential to match observations. However, it should be 
noted that the i?-charge chemical potential we are investigating is not the same as the more physically 
relevant chemical potential related to non-zero baryon number density. 
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We work in the framework of D = 5, A/" = 2 gauged supergravity, which is dual 
to A/" = 1 super- Yang Mills theory. In particular, we are interested in supersymmetric 
higher derivative terms, which have a highly constrained structure^. The four-derivative 
corrections to the leading order supergravity include a mixed gauge-gravitational Chern 
Simons term A ATt {R A R). The supersymmetric completion of this term was done in 
[19I . where an off-shell action was obtained for D = 5, Af = 2 gauged supergravity at 
the four-derivative level. In 20] we derived the corresponding on-shell Lagrangian, found 



corrected i?-charged black hole solutions, and studied their thermodynamic properties. 



We will use many of the results of 20|] to compute the shear viscosity. Our main result 
is that turning on /2-charge not only leads to violation of the bound, but enhances the 
effect, pushing rj/s further below 1/4tt. Furthermore, while the dependence of r] and s 
individually on the i?-charge is quite complicated, the ratio rj/s is remarkably simple. 

The general picture that emerges from such studies is that if we are interested in 
describing properties of the QGP (or other strongly coupled systems), we can try tuning 
the parameters available to us (whether N, A or the chemical potential), as long as we 
remain within the regime of validity of the supergravity approximation. Moreover, it is an 
interesting fundamental question whether violations of the bound can be related to any 
constraints on the dual gravitational side or consistency requirement of the underlying 
string theory. For instance, one may be able to relate the violation of the rj/s bound to 



the weak gravity conjecture of 21[, according to which there should be some states whose 
M/Q ratio is below the BPS bound. While this is an interesting avenue to explore^, the 
solutions that we have considered do not admit a nice extremal BPS black hole limit (since 
the extremal solution is the superstar geometry, with a naked singularity), and therefore 
do not lend themselves easily to such an analysis. 

The structure of this paper is as follows. In the following section, we present the 



^ Four derivative corrections in the presence of a chemical potential have been partially discussed in 
where and J"* corrections were considered, respectively. The supersymmetric Lagrangian, 
however, has RF^ and VFVF-type terms as well which were not previously considered. 
^ See [22| for investigating the effect of higher derivatives on the weak gravity conjecture. 
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on-shell four-derivative action and write down the non-extremal i?-charged background. 
Using this as a starting point, we then compute the shear viscosity in Section IIIII and 
conclude with a brief discussion in Section IIVI Some of the intermediate expressions are 
relegated to an appendix. While this work was being completed we became aware of 231], 
which overlaps with our results. 



II. M = 2 GAUGED SUPERGRAVITY AND i?-CHARGED BLACK HOLES 



Our starting point is five-dimensional M = 2 gauged supergravity. The physical fields 
in this theory are the metric g^j_y, graviphoton and gravitino ip^j^. The supersymmetric 
four-derivative corrections were obtained in [l^ using the superconformal tensor calculus 
methods worked out in 2J, |25|, l26|, l27|. By integrating out the auxiliary fields, the 
Lagrangian may be put into the form 20|] ^ 



+ - 



C2 

9' 



1 



16^3 



12^3 



16 



—RF^ + - F^,W,F^P + - V^F'PV^F,, + - V^F'^V^F, 



pp 



32\/3 



+—F^,F''PF,xF^^ 



64' 



41 
2304 



{F 



2\2 



(4) 



The four-derivative corrections are determined in terms of a single new dimensionless 
parameter C2 fcorrespondirig to coq^/24 in the notation of [2^). Holographic computation 



of the Weyl anomaly 13 



28 



29 



30| allows G5 and C2 to be expressed in terms of the 



anomaly coefficients a and c of the dual A/" = 1 gauge theory. This was worked out in 
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201 ]. with the result 



71 



g G5 = — . 

^ ^ 8a' 



c — a 



C2 



(5) 



We follow the conventions of 



19| and take [V^, V^]i;'^ = R^^p " yf and Rab = R, 



'ac b ' 
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Nonextremal i?-charged black hole solutions to the lowest order M = 2 gauged super- 



gravity were found in 3l|, and the corrections linear in C2 were worked out in 20|]. Using 
a parameterization convenient for the shear viscosity calculation, the flat-horizon black 
holes are given by the metric 

Hiu) 



ds 



2 9^rl \ f{u) 



and the gauge field 



At = gro] 



u VHiu) 



^3(1 + 



;dt - H{u)dx- 



-du 



(6) 



1 - 



l + qu 2^^ {1 + quY 



(7) 



The metric functions f{u) and H{u) are given by 

/ = (1 + quf - (1 + g) V + C2 

_ 9 

C2 



H 



u 



i + ,«-^,(i + ,)'^^^. 



(8) 



The above solution is fixed in terms of two parameters, tq (related to non-extremality) 
and dimensionless q (related to the i?-charge). At the two-derivative level, the horizon is 
located at u = 1, while the boundary of AdSs is at m = 0. At linear order in C2, however, 
the horizon location gets shifted to 

C2 (l + g)(3-26g-f 3g2) 



1 + 



12 



2-g 



(9) 



The temperature and entropy density were obtained in 20| 



T 



^Vo(2-g)(l + g) 
27r 

[grofH + qf/^ 



1/2 



1 - 



C2 10 - 59g - 4g2 - Sg^ 



4G', 



1 + 



{2-qf 
C2 21 + 14g-3g2" 
¥ 2-q _ 



(10) 



Note that, for g = 0, we may write the entropy density in terms of the temperature as 



2T(^a 



9 c — a 
1 + 4 — 



:ii) 



where we used the holographic relations ([5]). This reduces to the familiar s = tt'^N'^T^ /2 
3^ for A/" = 4 SYM, where a = c = N'^/A. 



III. COMPUTATION OF THE SHEAR VISCOSITY 



We compute the shear viscosity using the Kubo formula, following the methods devel- 
oped in In particular, we introduce a scalar channel perturbation to the metric 

9xy ^ 9xy ~l~ hxy, (12) 

where, for convenience, we define h^y = (j){t,u,x). Expanding the Lagrangian (jlj) to 
second order in the perturbation yields 



+Efy'L, + J , (13) 



where the fourier components of (f) are defined by 

0(t,u,x)= ! d^xdt<Pk{u)e'^^'^~^'\ (14) 



We note that this parameterization of the action with coefficients A, . . . ,F was originally 



used in j? I to handle the correction of IIB supergravity. However, it is general enough 
to accommodate the present case. The coefficients are even functions of the momentum, 
and are given explicitly in the appendix. 

Varying this action with respect to yields a fourth order differential equation. How- 
ever, since the higher derivative terms are multiplied by C2, we may reduce the order of 
the equation by working perturbatively in C2- To see this, we first consider the lowest 
order equation of motion 

where we have defined the dimensionless frequency 

,2 



0"+ f -- + = 0, (15) 



The lowest order metric functions 

/o= (l + gM)'-(l + g)V, Ho = l + qu, (17) 



are obtained by setting C2 = in ([8]). Taking additional derivatives of f[T5|) allows us to 
eliminate 0'" and 0"" terms in the full equation of motion. The result is rather simple: 

Notice that the form of this equation is almost identical to that of (fT5|) . the lowest order 
equation of motion, modified only by the presence of the corrected metric functions / and 
H as well as one new term, which is explicitly 0{c2)- 

Since the function f{u) vanishes linearly at the horizon m+, the point m = m+ is a 
regular singular point of the equation of motion f|T8|) . This suggests that we write 



0(«) = /(«rF(«), (19) 

where F{u) is assumed to be regular at the horizon. The exponent v is then obtained 
by solving the indicial equation. In the hydrodynamic limit, the lowest order solution is 
known Q, [l^ and is given by: 



0o = /o(«)''°<!l-y 



where 



Aln— — — —+3\n[l + {ai + l)u-a3U ) 

(^ + «! + 1 - 2q;3m)(^ - ai - 1) 



(20) 



ai = 3q, a2 = 3q\ = q\ E = {1 + q){l + AqY/\ A = -3^. (21) 



The exponent uq is given by 



(2-g)"(r+g)V2 



^0 = -TT^-XTT—Xm^ (22) 



and may be re-expressed as uq = —iu/A-KTo, where Tq is the lowest order temperature 
given in ffTOl) . Note that we have chosen incoming wave boundary conditions at the horizon 
as appropriate to the shear viscosity calculation. 

Adding higher derivative terms will have two effects on this solution, one being a 
correction to the function F{u) and the other a modification of the exponent u defined 
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above. For the exponent, solving the indicial equation gives 

C2 10 - 59g - 4g2 - 3q^ 



1 + 



47rT' 



(23) 



(2-g)(l + g)V2V 8 (g 

where the relation to the temperature ffTOj) is valid to linear order in C2- We may now 
substitute 0(m) = f{uYF{u) into the equation of motion (fTSl) and linearize in C2 to obtain 
an equation for F{u). While this is difficult to solve exactly, since we only need a solution 
in the hydrodynamic regime, it is sufficient to work to first order in uj (or equivalently v) . 
The solution for F{u) is quite complicated and can be found in the appendix. 

Given this solution, it remains to evaluate the on-shell value of the action. As explained 
in [3], the bulk action (fT3|) must be paired with an appropriate generalization of the 
Gibbons-Hawking term. In general, the fourth order equation of motion yields a boundary 
value problem for the two-point function where additional data must be specified {e.g. 
fields and their first derivatives at the endpoints). However, when working perturbatively 
in C2, the equation of motion reduces to a second order one, given by f|T8l) . This allows us 
to use a generalized Gibbons- Hawking term of the form 

F 



K, = —Acpk 



+ E{pi(f)'f. + 2po0fc)0- 



fc5 



(24) 



where 



Pi 



fo 



1 

u 



P2 



(25) 



are the coefficients in the lowest order equation of motion f|T5l) . 

Evaluating the on-shell action then amounts to evaluating a boundary term 



S 



(27r) 



7^k 



(26) 



where 



B-A 

+E<p 



El. 



l^'-k 



)0^0_, + i(C-A' 
/o u 



E'(f)l(f)_k 



E(p 



III A 



k + '^E 



ri2 3 



2 (f^k<P-k 



(27) 
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In order to compute the shear viscosity we need only the hmit of the above action as u 
approaches the AdS boundary {i.e. m ^ 0). It turns out that only the first and third 
terms contribute. This yields a value for the shear viscosity via the Kubo relation 



Tj = lim — lim(2ImjFfc) 

a;-^0 to M— ►O 



C2 5 + 6g + 



(2^ 



IGTrGs V ' 8 2-q 

Finally, dividing this by the entropy density flTI]]) gives a value for the shear viscosity to 
entropy density ratio of 



1 

s 



1 

47r 



l-C2(l + g) 



1 

Att L 



1 



c — a 



1 + 



(29) 



where we have rewritten C2 in terms of the anomaly coefficients c and a using 



IV. DISCUSSION 



The expression for t]/s, given in ( l29l) . is surprisingly simple, given that both t] and s are 
individually rather more complicated functions of the parameter q. This is presumably 
related to some form of universality, which holds even in an i?-charged background^. It is 
instructive to examine the contribution of the various terms in the four- derivative action 
to the result ( !29l) . We find that only four terms in (jll) are important. Writing 

1 



-R 



4 



+ 



C2 

9' 



we may arrive at the result 

7] 1 

s Arc 

Note that setting to their actual values in (jlj) reproduces (!29|) . 



1 — 4c2(2ai — q{ai + 60:2 — 60:3 + 3a4)) 



(30) 



(31) 



^ Of course, the simplest result possible would have been to obtain tj/s independent of q. But this is 
clearly not the case here. 
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The shear viscosity to entropy density ratio was independently derived in 23|, where 
it was found to depend only on terms explicitly involving the Riemann tensor [i. e. the ai 
and a2 terms in ( l30i) ]. This appears to differ from the result found above. However, by 
the use of Bianchi identities and integration by parts we can cast the gradient terms into 
the form 



(2a3 - 04) 



- F^.V^V.F'^" + F'^pF.^R: - \Ra,,„F>'' FP^ 



(32) 



The first two terms do not contribute to the rj/s ratio, while the last term will add to the 
original 0:2 term to give an effective a2 = 02 — 0^3 + 0^4/2, so that (13T1) may be rewritten 
as 



1 

s 



1 



(33) 



1 - 4c2(2ai - q{ai + 60:2)) 

This agrees with the result of [231] provided the difference in signature conventions is taken 
into account. 

Finally, we return to the M = 1 SYM shear viscosity result of fl29p . In order to express 
this in terms of physical quantities, we wish to relate the parameter q to the i?-charge 
chemical potential and temperature. Since q only enters into (!29|) at the next-leading 
order, we can use the leading order expressions in pinning down q. The chemical potential 
for i?-charge $ is identified as the difference of At between horizon and boundary 33l. l34l|. 
At lowest order, ([7]) yields 

$ = gro^?,q{l + q). (34) 



Comparing this to the temperature 

To = 

allows us to write 



9^ 
27r 



(2-g)(l + g) 



1/2 



2$2 



1 + -<i>^ - 




(35) 



(36) 
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where $ = g^/27TT is the dimensionless chemical potential. Note that q is an increasing 
function with respect to $, with q = when $ = 0. The possible value of q ranges as 

< g < 2. (37) 

Substituting (l36l) into ( l29i) then gives 

9. ( 4 _ _ \ 8 _ \ \ 

(38) 



?7 1 

s Atc 




2$2 

Since q is non-negative, this demonstrates that turning on an i?-charge chemical potential 
only increases violation of the t]/s bound, provided c — a > 0. Taking the range (1371) into 
account, we see that adjusting the i?-charge yields a range of values 

where we have again assumed c — a > 0. 

In conclusion, we have explored the effect of a background i?-charge on the shear 
viscosity to entropy density ratio rj/s. While the leading order ratio rj/s = I/Att is 
universal, i?-charge corrections do turn up at the order. For known theories with 
a holographic dual, where c — a > 0, the conjectured 1/Att bound is generally violated 
for arbitrary chemical potential. We caution, however, that this is a parametrically small 
violation appearing at 0{1/N) in the large N limit. In principle, it would be desirable to 
obtain a more robust result. However, this is hindered by difficulties in obtaining exact 
solutions to the full equations of motion {i.e. beyond the linearized limit). While this 
can be done in certain cases such as Gauss-Bonnet gravity, the natural supersymmetric 
organization of the higher derivative Lagrangian (jl]) is not of this form. It would be 
interesting to see if a modified universality relation for rj/s can be obtained for arbitrary 
forms of the higher derivative gravity theory. 
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VI. APPENDIX 

The quadratic action for the scalar channel perturbation is given in f|T3l) in terms of 
six coefficients A, . . . ,F. Here we present their explicit forms: 



Aiu) 
B{u) 



C{u) 



D{u) 



E{u) 
F{u) 



-fo + C2 

u 



uj^Hl 2ufQ{l + qf{hqu-l) 32^ V (1 + g^u^{l + qf 

7^ ^ M 3 ^ Ho 



3/o 



u 



+ C2 



^%r2 g^iMu + lfHl ^(^^(l + g)3(56gV + Tgii + 11) 



3^ 



6 



^V(l + g)6(26gV- 17gM+17) „ 2,. ^ ^3 2^ 3^7^(1 + 

8(7 (1 + q) qu + - 



2g'^{Aqu - ?,)Hl 2g\l + qf{2qu + 1) 



Hn 



+C2 



£^ ((4gM + l)i/o' - (1 + g)'(-llgw=^ + Uu'')Ho) 



6m/( 

5(2(1 + g)3(4gV + 45^^^ + 3) g^u^{l + g)6(4g3^3 _ 7g2^2 _ ^2qu + 15) 



3//0 

2g^Hl - g^qu^{l + g)^ 



2H^, 



2 -^0 



+ C2 



a; 



^ + ""'^'is ft (2(3l9« - 9)//o' - 3n^(l + g)3(5gV - Aqu + 11)) 



12< ' 48/5 
19^2^(1 + qf 3^2^(1 + qf{Qq^u'^ - llqu + 1) 



ml 



2HI 



C2 



C2/0 



2(2(4gM + l)Hl - u^{l + qfilqu + 4)) 



ml 



(40) 



Here we also present the 0{c2) solution for 0. Writing 0('u) = f{uYF{u), we may 
expand F{u) to first order in both C2 and u; 



F{u) = Fo{u,uj) + C2{Fw{u) + LuFuiu)). 



(41) 



13 



Since F{u) satisfies a second order equation (after linearizing in C2 and using tlie lowest 
order equation of motion), it is consistent to clioose tlie boundary conditions sucli tliat 
F{u) is normalized at the boundary {F{0) = 1) and is regular at the horizon. 

The function Fq{u,uj) is given by the expression in the curly brackets in (1201) . while 
the remaining functions are 



Fioiu) 
Fniu) 



0, 

(1 + g)3/2(llg5 + 4^4 ^ ^79^3 _ iQ^2 _ _ 15) 
32g2(l+g)2(g_2)3 

i{q + l)3/2(60g6 + 99^^ + 648g^ - 69q^ - 15V - 104g - 16) 



i \n{q W — 3qu — u — 1) + 



TV 



16(4g + l)3/2(g + l)2(g-2) 



tanh ^ 



:i + 3g) 



tanh 



_i 2q^u-{l + 3q) 



(4g+ l)i/2(g + i) (4g+ 1)1/2(^ + 1^ 

i ln(l + gM)(l + g)3/2 i(g + l)3/2(_4g5 ^ 21g4 + 143^3 _ 21^2 _ 39^ _ 



8g2 8g4(4g+ l)(g-2)2 

i(g + l)3/2(4g7 _ 27g6 + 64^^ + 511g^ + ISTg^ - 128g2 - 57g - 6)9^2 

8(1 + gM)g4(g3y2 _ _ ^ _ 1)^4^ + l)(g - 2)2 
i{-12q^ + 102g5 + 605g^ + 63q^ - YJlq^ - 63g - 6)m 



8(1 + qu)q\q^u^ 



3qu-u- l)(4g + l)(g - 
i(4g5 + 21g^ + U3q^ - 21q^ - 39q - 6) 
8(1 + qu)q^{q^u^ - 3qu - u - l)(4g + l)(g - 2)^ 



(42) 
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